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ABSTRACT 

Chu,  P.C.,  Garwood,  R.  and  Muller,  P.,  1990.  Unstable  and  damped  modes  in  coupled  ocean  mixed  layer  and  cloud  models.  J. 
Mar.  Syst.,  1:  1-11. 

The  ocean  mixed  layer  and  clouds  are  coupled  by  the  fluxes  of  momentum,  heat,  and  water  mass  at  the  interface.  The 
importance  of  the  fluxes  of  momentum  and  heat  is  well  recognized  by  both  meteorologists  and  oceanographers.  However,  the 
water  mass  flux  has  been  given  considerable  attention  only  in  atmospheric  models  since  the  latent  heat  release  is  an  important 
source  of  energy  for  the  atmospheric  general  circulation.  The  water  mass  flux  is  given  less  attention  in  ocean  models  although 
it  is  realized  that  evaporation  and  precipitation  contribute  to  the  surface  buoyancy  flux  which  influences  the  depth  of  mixing 
and  the  thermohaline  circulation. 

Clouds  and  the  ocean  mixed  layer  are  coupled  by  both  the  heat  and  moisture  fluxes  across  the  air-ocean  interface.  Two 
time  scales  are  demonstrated  in  this  paper:  a  sea  surface  temperature  (SST)  evolution  time  scale,  tt.  that  is  virtually  controlled 
by  the  oceanic  planetary  boundary  layer  (OPBL)  and  a  cloud-SST  coupling  time  scale,  t„  t.  These  two  time  scales  depend  on 
the  stability  of  the  marine  atmospheric  boundary  layer  (MABL).  The  more  unstable  the  atmosphere,  the  shorter  the  time 
scales  will  be.  For  a  stable  atmosphere,  iy-1-3  years,  and  r„  T~  0.3-1  years.  However  for  an  unstable  atmosphere, 
Tr  -  20-30  days,  and  r„  r  —  3-6  days. 

An  air-ocean  coupled  model  is  presented  in  this  paper  for  two  different  regimes:  (1)  the  non-entraining  ocean  mixed  layer 
case  and  (2)  the  entraining  mixed  layer  case.  The  model  results  demonstrate  that  the  exchanges  of  heal  and  water  across  the 
sea  surface  lead  to  both  growing  and  decaying  modes  of  oscillation  on  the  two  time  scales  due  to  the  stability  of  the 
atmosphere.  These  oscillatory  solutions  are  entirely  thermodynamic  and  do  not  require  wave  dynamics  for  their  existence. 


Introduction 

Since  1970  significant  progress  has  been  made 
both  in  our  ability  to  carry  out  air-sea  interaction 
field  work  and  in  our  understanding  of  many  of 
the  processes  found  on  both  sides  of  the  air-sea 
interface.  Many  studies,  both  observational  and 
theoretical,  have  shown  that  the  surface  wind  and 
the  SST  are  two  important  elements  in  the  air-sea 
coupled  system  (Chu,  1989). 

The  feedback  mechanism  between  clouds  and 
ocean  mixed  layer  can  be  explained  as  follows. 
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First,  clouds  reduce  the  incoming  solar  radiation 
at  the  ocean  surface  by  scattering  and  absorption, 
which  cools  (relatively)  the  ocean  surface  layer  by 
increasing  mixed  layer  entrainment.  The  cooling 
of  the  ocean  mixed  layer  lowers  the  evaporation 
rate,  which  will  diminish  the  clouds.  This  is  a 
negative  feedback  mechanism.  Second,  precipita¬ 
tion  dilutes  the  surface  salinity,  stabilizing  the 
upper  ocean  and  reducing  mixed  layer  deepening. 
The  mixed  layer  may  be  caused  to  shallow  if  the 
downward  surface  buoyancy  flux  is  sufficiently 
enhanced  by  the  precipitation.  The  reduction  in 
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Fig.  1.  Main  physical  processes  in  the  two  adjacent  boundary  layers. 


mixed  layer  depth  will  increase  the  sea  surface 
temperature  by  concentrating  the  net  radiation 
plus  heat  fluxed  downward  across  the  sea  surface 
into  a  thinner  layer.  The  increase  of  SST  augments 
the  surface  evaporation,  which  in  turn  produces 
more  clouds.  This  is  a  positive  feedback  mecha¬ 
nism.  Figure  1  shows  the  main  physical  processes 
and  fluxes  of  heat,  mass  and  momentum  in  the 
OPBL  and  the  marine  atmospheric  boundary  layer. 
Figure  2  illustrates  the  feedback  pathways  be¬ 
tween  the  clouds  and  the  ocean  mixed  layer  (Chu 
and  Garwood,  1989). 

The  feedback  mechanism  between  clouds  and 
ocean  mixed  layer  largely  depends  on  the  dy¬ 
namics  of  ocean  mixed  layer.  Ocean  surface 
buoyancy  flux  either  generates  or  dampens  the 
turbulent  kinetic  energy  in  the  upper  ocean  de¬ 
pending  upon  the  direction  of  the  flux.  Upward 


Fig.  2.  Feedback  paths  between  OPBL  and  clouds. 


buoyancy  flux  tends  to  cause  the  water  column  to 
be  statically  unstable,  which  generates  turbulence. 
On  the  otherhand,  downward  buoyancy  flux  makes 
the  water  column  to  be  statically  stable,  damping 
the  turbulence.  In  the  case  of  upward  buoyancy 
flux,  or  in  the  case  when  the  shear  production 
exceeds  the  buoyant  damping,  the  mixed  layer 
entrains  water  from  below,  increasing  the  poten¬ 
tial  energy.  If  there  is  sufficiently  strong  down¬ 
ward  buoyancy  flux,  the  shear  production  may  be 
insufficient  to  overcome  the  buoyant  damping  and 
turbulence  cannot  reach  the  former  mixed  layer 
depth,  causing  mixed  layer  retreat  (damping  reg¬ 
ime).  It  is  expected  that  the  effect  of  clouds  on  the 
ocean  mixed  layer  is  different  for  these  two  oc¬ 
eanic  regimes. 

Since  clouds  have  significant  effects  on  the 
large-scale  atmospheric  circulation  through  the 
transfer  of  heat,  moisture  and  momentum  and  on 
the  ocean  mixed  layer  through  the  attenuation  of 
solar  radiation  at  the  ocean  surface,  and  since  the 
SST  is  an  important  factor  for  the  development  of 
clouds,  the  feedback  mechanism  mentioned  above 
has  a  potentially  significant  importance  for  air-sea 
interaction,  weather  and  ocean  prediction. 

Although  our  coupled  model  is  one-dimen¬ 
sional,  we  are  aware  of  the  importance  of  horizon¬ 
tal  advection  and  the  limitation  of  one-dimen¬ 
sional  models.  However,  the  intent  of  this  work  is 
to  develop  a  formalism  to  examine  thermody¬ 
namic  feedback  between  the  two  fluids.  Because 
we  wish  to  concentrate  on  the  thermodynamic 
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interaction,  horizontal  advection  is  ignored  ini¬ 
tially. 


Atmospheric  model 


Atmospheric  surface  layer 


Using  similarity  theory,  if  a  variable  (wind, 
temperature,  or  moisture,  etc.)  is  appropriately 
scaled,  then  its  vertical  profile  follows  a  universal 
function  whose  form,  in  general,  may  be  de¬ 
termined  empirically.  For  the  barotropic  atmo¬ 
sphere  the  appropriate  scales  for  wind  speed  and 
boundary  layer  height  are 

U=Ut,ha~jf-\  (1) 


where  Ug  is  the  geostrophic  wind  speed,  which  is 
assumed  here  to  be  10  m/s,  /  is  the  Coriolis 
parameter  and  wa„  is  the  atmospheric  friction 
velocity.  Using  the  similarity  functions  a,  b  and  c, 
the  geostrophic  drag  coefficient  Cg  is  computed  by 
Yamada  (1976)  as 


nu. 

I  £41 


(2) 


where  k  is  the  von  Karman  constant,  0.4.  The 
roughness  parameter,  z0,  is  approximately  1.5  X 
10  4m.  The  ratio  /ia/z0  is  0.6  X  107  for  an  MABL 
height,  /ia  =  1  km.  The  heat  and  moisture  transfer 
coefficients  are  determined  by 


Ch  — 


w'0' 


«a*[^~0Uo)] 


,  CE  — 


Wa'q'lo 


“*•[$- q(zo)] 

(3) 


where  w,  8 ,  q  are  the  vertical  velocity,  potential 
temperature  and  specific  humidity,  respectively. 
Here  the  symbol  “A  ”  denotes  the  values  being 
taken  at  the  top  of  the  MABL.  The  heat  transfer 
coefficient  is  given  by 

r  =  JL 
Ch  Pr0 

where  Pr0  is  the  turbulent  Prandtl  number  for 
neutral  stability,  having  a  value  of  0.74  according 
to  Businger  et  al.  (1971).  Based  on  the  Wangara 
data,  the  similarity  functions  a,  h,  c  are  experi¬ 


mentally  determined  for  the  stable  atmosphere 
(Yamada,  1976): 

3  =  1.855  -  0.38/ia/La  (0  «  hJL&  <  35), 
a  =  —2.94 (h,/L,  -  19.94)1/2  (35  <  /r,/La); 
b  =  3.02  +  0.3 hJL,  (0  <  hJL.A  35), 
b  =  2.85(  hJL,  -  12.47)1/2  (35  <  hJLa ) ; 
c=  3.665  -  0.829/ia/La  (0  «  h.JLa  <  18). 
c=  — 4.23(/ia/La  —  1 1  -21  )1/2  (18  <  /ta/La)  (5a) 
and  for  the  unstable  atmosphere  ( hJL a  <  0): 

10.0-  8.145(1  - 0. 008376// yL,)-173, 
b  =  3.02(1  -  3.29 /ia/La)_1/3, 
c=  12.0-8.335(1  -  0. 031 06/i a/La)_,/3  (5b) 

where  La  is  the  atmospheric  Obukhov  length  scale. 
It  is  noteworthy  that  the  formulae  (5a,b)  were 
obtained  form  the  Wangara  atmospheric  boundary 
layer  over  land  rather  than  from  the  MABL.  There 
have  been  no  comparable  observations  in  the 
MABL  as  yet.  Nevertheless,  we  assume  that  these 
similarity  functions  (5a,b)  are  valid  for  the  MABL. 
Furthermore,  since  the  computation  method  for 
the  moisture  transfer  coefficient  is  not  well  estab¬ 
lished,  in  this  paper  we  assume  that 

Ce  =  Ch  (6) 

Substitution  of  (5a,b)  into  (2)  and  (4)  leads  to  the 
apparently  strong  dependence  of  Cg,  CH,  CE  on 
the  atmospheric  stability  parameter  h.JLa,  as 
shown  in  Fig.  3.  These  parameters  have  much 
larger  values  for  the  unstable  atmospheric  than  for 
the  stable  atmosphere,  i.e., 

Cg~  0.0316,  CH,  CE~  0.063  (ha/La<0) 

Cg  ~  0.003,  CH,  CE  -  0.004  {h.JLa  >  0)  (7) 

Simple  cloud  model 

The  time  rate  of  change  of  cloud  cover  is  pro¬ 
portional  to  the  moisture  supply  divided  by  the 
amount  of  water  vapor  necessary  to  produce  the 
model  cloud.  The  main  processes  causing  the  cloud 
dissolution  are  precipitation  and  mixing  with  the 
environmental  air.  The  cloud  evaporation  due  to 
mixing  with  ambient  air  is  a  complicated  problem. 
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Fig.  3.  Dependence  of  Cg,  CH.  CE  on  the  atmospheric  stability 
parameter  ( h  a  /La ). 

and  is  neglected  for  the  sake  of  simplicity  here. 
Thus  the  equation  for  cloud  cover  is  reduced  to 
9 n  _  (Ml0  +  E-  Pr) 

dt  hc  W 

where  hc  is  the  total  amount  of  water  vapor 
needed  to  create  the  cloud  over  a  unit  area,  Ml0  is 
the  large-scale  horizontal  moisture  convergence  in 
the  column  of  atmosphere  per  unit  area,  E  is  the 
surface  evaporation  rate  and  P  is  the  precipitation 
rate.  From  mean  distributions  of  temperature  and 
mixing  ratio  in  the  environmental  air  outside  the 
cloud  and  inside  a  deep  cumulus  cloud  (Kuo, 
1965),  we  estimate  that  hc  —  5  cm. 

Relationship  between  precipitation  rate  and  cloud 
cover 

By  linear  regression  of  hourly  rain  amounts  and 
satellite  IR  brightness  data  obtained  during  Phases 
I,  11  and  III  of  GATE.  Albright  et  al.  (1985) 
suggested  a  linear  relationship  between  average 
precipitation  rate  Px  in  boxes  on  a  side  and  cloud 
cover  n  of  the  boxes  by  clouds  with  tops  colder 
than  -36°C: 

/>r(ms*')  =  (0.472  +  8.333n)XlO"7  (9) 

This  result  verified  Arkin's  (1979)  earlier  analysis 
for  the  GATE  B-scale  array. 

Cloud  effects  on  the  net  radiation  at  the  ocean 
surface 

Clouds  reduce  the  solar  radiation  incident  at 
the  ocean  surface  by  scattering  and  absorption. 


These  effects  are  included  in  Budyko’s  (1978)  for¬ 
mula 

Rs=[l-asnn-aso(l-n)]Rs0  (10) 

Here  Rs0  (340  W  m  2)  is  the  solar  radiation 
absorbed  by  the  ocean  surface  layer  under  a  clear 
sky.  The  parameters  asn  and  as0  represent  albedos 
of  the  earth-atmosphere  system  with  complete 
cloud  cover  and  a  cloudless  sky,  respectively,  and 
have  the  following  values: 

«Sn  =  0-46,  «so  =  0.2 

The  ocean  surface  emits  longwave  radiation  to  the 
atmosphere  and  to  space.  However,  clouds,  as  well 
as  dry  air,  partially  absorb  the  radiation  and  re¬ 
emit  longwave  radiation  back  to  the  ocean  surface. 
Thus  the  net  upward  energy  loss  by  longwave 
radiation  at  the  ocean  surface,  Rh ,  is  corrected  for 
the  downward  radiation  by  the  clouds  and  the  air. 
From  longwave  radiation  data,  Budyko  (1978)  de¬ 
rived  a  semi-emprical  formula: 

Rb  =  a  +  bT- (o,  +  bx)n  (11) 

The  dimensional  coefficients  a,  h ,  a,  and  6,  are 
a=  -377.6  Wm~2,  h  =  2.2  Wm“2K“', 
at  =  -389.8  Wm~2,  b ,  =  1.6  Wm'2K_1 

The  effects  of  clouds  on  buoyancy  flux  at  ocean 
surface 

For  the  case  of  an  ocean  surface  without  ice 
(low  and  middle  latitudes)  the  surface  buoyancy 
flux,  B,  has  two  components: 

B=--^~-Pg(E-Pr)S  (12) 

rwl  pw 

Here  a  is  the  sea  water  thermal  expansion  coeffi¬ 
cient,  fi  is  the  salinity  contraction  coefficient,  pw 
is  the  water  density  and  cpw  is  the  sea  water 
specific  heat  under  constant  pressure.  The  surface 
heat  flux,  F  (upward  positive),  is  computed  by 

F=Rh-  Rs  +  LpwE  +  Hs  (13) 

where  L  is  the  latent  heat  of  vaporization  and  Hs 
is  the  sensible  heat  flux  to  the  air.  The  variables 
T_h  and  S_h  are  the  temperature  and  salinity  of 
the  water  immediately  below  the  mixed  layer  that 
may  be  entrained  into  the  mixed  layer.  Standard 
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bulk  formulae  are  used  to  calculate  the  surface 
evaporation: 

E  =  paCDt/gCE  [ q% (Tw)  -  q] /pw  ( 14) 

and  the  sensible  heat  flux  from  the  ocean  surface: 

f/s  =  pacpacDf/gcH(rw-fl)  (15) 

where  qs(T)  is  the  saturated  mixing  ratio,  cpa  is 
specific  heat  of  the  atmosphere,  and  6  and  q  are 
the  air  temperature  and  mixing  ratio  at  the  top  of 
the  MABL. 

The  effect  of  clouds  on  the  buoyancy  flux  at 
the  ocean  surface  is  two-fold:  (1)  decreasing  B 
through  the  increase  in  the  net  heat  loss  at  the 
ocean  surface.  F,  by  reducing  the  incoming  solar 
radiation  and  (2)  increasing  B  due  to  precipita¬ 
tion. 


Ocean  mixed  layer  model 


Most  models  which  include  thermodynamic  ef¬ 
fects  regard  the  upper  layer  as  a  well-mixed  turbu¬ 
lent  boundary  layer  which  exchanges  heat  and 
moisture  with  the  atmosphere  and  entrains  water 
from  below.  The  heat  and  salinity  equations  take 
the  forms 


9  T  F 

=  -we(rw-r_h) - ^+at  (16) 

Awf  -  -wc(S-S_h)  +  (£-Pr)S  +  /ls  (17) 


where  hw  is  the  mixed  layer  depth  and  AT  and  As 
are  the  horizontal  advection  for  temperature  and 
salinity,  respectively.  The  entrainment  velocity  is 
wc  and  is  parameterized  as  (Chu  and  Garwood, 
1988): 


=  (C,u3w.-C2Bh) 

We  J  gf>[a(Tlt-r_/l)-/i(S-S_,)] 


(18) 


where  C,  and  C2  are  tuning  coefficients  and  «w» 
is  the  water  surface  friction  velocity,  which  is 
computed  by 


for  an  atmospheric  friction  velocity  of  30  cm/s, 
mw*  -  1.1  cm/s.  The  symbol  A  is  a  Heaviside 
function  of  (C:u^.  —  C2Bh).  When  (C,m*.  — 


C2Bh)  >  0  there  is  sufficient  turbulent  kinetic  en¬ 
ergy  to  entrain  and  mix  water  from  below  and 

A  =  1 

which  represents  the  entrainment  regime.  The  en¬ 
trainment  rate  is  determined  by  (18)  and  sub¬ 
stituted  into  the  following  equation  to  prognosti¬ 
cate  the  mixed  layer  depth  h. 

dh 

_  =  we-w_A  (20a) 


Here  w_A  is  the  vertical  velocity  at  the  mixed 
layer  base. 

When  (C,w*  —  C2Bh)  <  0  there  is  not  enough 
turbulent  kinetic  energy  to  entrain  water  from 
below,  wc.  is  set  to  be  zero,  i.e., 

A  =  0 


which  is  called  the  surface  damping  regime.  The 
mixed  layer  depth  is  calculated  diagnostically  from 
a  balance  of  the  remaining  terms  in  (18)  and  it 
equals  the  oceanic  Obukhov  length  scale,  Lw: 


C,mw. 

C2B 


(20b) 


One  example  for  the  surface  damping  regime  is 
the  western  Pacific  warm  pool  region  where  the 
excess  precipitation  effect  prevails  over  the 
buoyancy  loss  effect  of  heat  loss  at  the  ocean 
surface,  causing  the  net  buoyant  flux  to  be  down¬ 
ward. 

The  temperature  and  salinity  equations  (16) 
and  (17)  describe  the  balance  of  storage,  entrain¬ 
ment  and  heating/water  mass  flux. 


Basic  equations  for  perturbations 

The  mean  state  of  the  coupled  cloud-OPBL 
system  (7/.,  5,  h ,  E,  Pr,  we)  is  evaluated  from  the 
steady-state  solutions  for  the  prognostic  equations 
(8),  (16),  (17)  and  (20a). 

When  the  coupled  system  is  perturbed  from  its 
equilibrium  state,  the  thermodynamic  feedback 
mechanism  between  the  cumulus  clouds  and  the 
oceanic  mixed  layer  causes  the  perturbation  to 
either  grow  (positive  feedback)  or  dampen  (nega¬ 
tive  feedback).  The  principal  purpose  here  is  to 
study  possible  one-dimensional  thermodynamic 
feedback  mechanisms  between  clouds  and  the  oc¬ 
eanic  mixed  layer,  within  the  limitations  of  the 
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simplifying  assumptions.  Hence  the  energy  ex¬ 
change  at  the  air-ocean  interface  is  a  primary 
focal  point.  Therefore,  we  shall  neglect  initially 
the  perturbations  of  those  variables  which  are  not 
directly  related  to  the  exchange  at  the  air-ocean 
interface,  such  as  the  perturbations  of  the  horizon¬ 
tal  advection  A'r  and  A's,  the  perturbations  of  the 
atmospheric  variables  at  the  top  of  the  MABL,  q’ 
and  O'  and  the  perturbations  of  the  oceanic  varia¬ 
bles  at  the  base  of  the  OPBL.  wlA,  T'_h  and  S'h. 
From  the  basic  equations  of  the  coupled  system 
(8),  (16).  (17)  and  (20a  or  b),  the  perturbation  of 
cloud  cover.  n\  is  given  by 

9«' 

“9 T 


1  I  dE  ,  9 P,  , 

77  -  — —n 


h  I977.‘w  0/i 

The  equations  for  77',  5',  are 

977  _ _ 1 

9/ 


(21) 


9 F  .8  F  ,  F 
HAJU8""  3r-  "  ’ 


w  pw 

-_A(xT-  +  f< 

'-wJ1rbs‘+bE'-p:) 

,  (Er~E)S  ,,  A  (W  ,  AS 

+ - p - h^-A^jS  +-p we 


(22) 


(23) 


The  equation  for  h'w  has  different  forms  for  the 
two  regimes.  For  the  entrainment  regime  it  has  the 
prognostic  form 
9  h' 

9/ 

where 

S3 


=  w. 


(24a) 


we'  A  b  =  -gw_h(aTs'  -  PS!)-  C,^fh'  +  C2£' 

(25) 

For  the  surface  damping  regime,  it  has  the  di¬ 
agnostic  form 


From  (9)  and  (14)  £'  and  P’  are  computed  by 


d  E 


£'  =  — —  T'  P’  = 

ftp  '  w  ’  1  r 


9^ 

3  n 


(26) 


Three  basic  time  scales 

Three  time  scales  are  found  from  (21 )— (24) 

(a)  Cloud  lime  scale 

J_9 Pj_ 

7 "  /ic  9/? 

Using  (9)  and  taking  hc  =  5  cm.  we  have 
r„  —  0.6  day 

(b)  SST  variation  time  scale 

1  9  £ 


(27) 


T r  '  = 


p  e  h  9  77. 

rwc  pw  n  w  w 


(28) 


Based  on  the  assumption  that  0  and  q  are  de¬ 
termined  by  the  large-scale  atmospheric  motion 
only,  from  the  heat  balance  (13)  we  have 

1L-  r'L  . 

9 Tw  +  r  RT}  CH  C^  H 


pa 


^2qs  CE 


Therefore 

,  =  ^  +  ^  , 
PwCpw^w  \  Cpa^yT’w  CH  / 

Cloud  and  SST  coupling  time  scale 

1  9£  9£  / 97s, 

Vr_  PwCpw^w  9,1  97W  9” 
Utilization  of  (27)  leads  to 


(29) 


B 

9  £  ‘ 

I  ag  9  £ 

4-  ft  c  9£  W' 

T~1  = 

B 

\  PwC  pw  9£w 

9£w  j 

rn.T 

Pw^pw  ^  W  h  c  ^  V 

-j8gS(Pr-£)S' 


«g  dF  -bPr  \  , 

+  W^'/!sS^1 


(24b)  9£  = 


CJOh 


Using  (10),  (11)  and  (13)  we  have 
9£  9 (*„-*.) 


(30) 


9n 


9n 


15Wm 
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Fig.  4.  Dependence  of  rT  and  t„  T  on  the  atmospheric  stability 
parameter  (/ia/Z.a). 


Among  these  three  times  scales,  t„  is  the  shor¬ 
test.  The  other  two,  tt  and  t„  r,  largely  depend  on 
the  parameters  Cg,  CH  and  CE,  which  are  func¬ 
tions  of  the  atmospheric  stability.  Figure  4  shows 
the  dependence  of  tt  and  r„  r  on  the  atmospheric 
stability.  For  the  unstable  atmosphere  (/ia/La  <  0), 
t„  r~  3-6  days  and  tt~  20-30  days.  For  the  sta¬ 
ble  atmosphere  (/7a/£a  >  0),  t„  7  -  0.3-1  years  and 
tt~  1-3  years. 

Since  the  time  scale  for  cloud  feedback  is  so 
much  shorter  than  that  for  SST  feedback,  i.e., 
t„  «  tt,  the  cloud  cover  perturbation,  almost 
instantaneously  follows  the  SST  for  the  tempera¬ 
ture  feedback: 


/ 

n 


ZE/dT„ 

a />r/a«  w 


(31) 


Equation  (9)  shows  that 

3  p 

-  8.33  X  10“7ms“\  (32) 

and  for  Tw  =  25  °C,  ua.  =  0.1  m/s. 


-  Cem3 .  d  VJ--  ^  ~  4.63  X  10~7ms  'K  1 


37; 

9£/97; 

9Pr/9;i 


drw 

-  0.56K-1 


(33) 


Therefore,  a  1°C  change  in  SST  implies  a  0.56 
change  in  cloud  cover. 


Surface  (lamping  regime 


Neglecting  the  small  terms  in  the  prognostic 
equation  for  7/  and  S’  and  eliminating  three 
among  the  four  variables  n' ,  Tw',  S'  and  h'w  from 
(22),  (23),  (24b),  (31),  we  obtain  the  second  order 
differential  equation. 


9^ 
3  r 


Y 


y-M 

MV2 


(y-m)2 


'  at 

)2*-o 


(34) 


where  ^  represents  7/,  5'  and  //v.  There  are 
two  nondimensional  parameters, 

/  PMcpw 

M  ^(Vr+rf1)  ’ 

4fe(Pr-£)5 
Y  Pghv,.S(r~^  +  Tf1) 

that  indicate  the  relative  importance  of  mean  heat 
and  mean  salinity  fluxes  in  the  mean  surface 
buoyancy  flux  B. 


Entrainment  regime 


In  the  entrainment  regime,  the  four  variables 
n\  7/',  S',  vv/  and  h'w  are  obtained  from  three 
prognostic  equations  (22).  (23)  and  (24a),  and  two 
diagnostic  equations  (31)  and  (25).  Eliminating 
four  variables  from  these  five  equations  we  obtain 
third-order  differential  equation: 

/a3  a2  a  \ 

(a?  +  flla7I  +  "237  +  a3r  =  0  (36) 


where  a 5,  a2,  o2  are  functions  of  the  following 
two  parameters  (Chu  and  Garwood,  1989): 


1  a  B 
~  ghc  3 n 


(37a) 


which  is  the  dependence  of  the  surface  buoyancy 
flux  change  on  cloud  cover  and 


_  ag  AT 
=  Ah 


(37b) 


which  is  the  fraction  of  the  mixed  layer  base 
density  jump  induced  by  temperature  to  the  total 
density  jump.  Here 

Ab  =  g(a  AT- P  AS) 


(38) 
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is  the  reduced  gravity,  and  A T  and  AS  are  the 
mean  temperature  and  salinity  jumps  at  the  base 
of  the  ocean  mixed  layer: 

Ar=  Fw-r_*,  as  =  s-s_a  (39) 


Solutions 


It  is  noteworthy  that  the  entrainment  regime 
(36)  is  of  the  third  order,  whereas  the  damping 
regime  (34)  is  of  the  second  order  because  of  the 
change  between  prognostic  and  diagnostic  form 
for  the  respect  of  mixed  layer  depth  equations 
(24a,b).  The  general  solutions  of  the  two  regimes 
have  the  form: 


'P='Ld,  exp  fat) 


(40) 


where  dt  {j  =  1,  2,  ...)  are  the  integral  constants 
and  Oj  (j  =  1,  2,  ...)  are  eigenvalues.  In  the  solu¬ 
tion  (40),  there  are  two  terms  for  the  surface 
damping  regime  and  three  terms  for  the  entrain¬ 
ment  regime.  For  the  surface  damping  regime,  the 
eigenvalues  are  the  roots  of  the  second-order  alge¬ 
braic  equation 


cr~  — 


+ 


+  r f1). 


f*Y" 


(y-m) 


;(T«.r  +  T7-1)2  =  0 


(41) 


For  the  entrainment  regime,  the  eigenvalues  are 
the  roots  of  the  following  third-order  algebraic 
equation: 

a'1  +  a,a2  +  a2o  +  =  0  (42) 


Instability  and  oscillation  criteria 


The  instability  criteria  for  the  thermodynami¬ 
cally  coupled  air-ocean  system  are 


I<  0  decaying 

=  0  neutral.  (43) 

>  0  growing 

where  a  is  a,.  a2,  or  a3,  the  roots  of  the  algebraic 
equations  (41)  or  (42).  The  oscillation  criteria  for 
the  coupled  system  are 


I  =  0  nonoscillatory 
Im(a)(  „ 

(  =£  0  oscillatory 


(44) 


Surface  damping  regime 


The  roots  o,  and  o2  are 
...  (2/*  —  y)  /  -i  .  -n 

°‘-2_  2(y  —  p)  '  T"  t  +  Tr  > 


MY' 


(2m -Yr 


(45) 


The  condition  for  the  generation  of  growing/de¬ 
caying  modes,  which  can  be  deduced  from  (43).  is: 


2^  ~  Y  I 
Y  —  M  I 


<  0  decaying 
=  0  neutral, 
>  0  growing 


(46) 


and  the  condition  for  oscillatory /nonoscillatory 
modes  is: 


/  2/t  \2(  >  p  nonoscillatory 

\  Y  /  \  <  p  oscillatory 

Separation  of  different  modes  in  the  y  —  p  plane 
is  shown  in  Fig.  5.  For  the  buoyant  damping 
regime  ( B  >  0,  i.e.,  y  >  p  >  0)  which  corresponds 
to  the  western  Pacific  warm  pool  regime  with 
strong  precipitation,  the  relative  importance  of  the 
mean  surface  heat  and  salinity  fluxes  in  the  mean 
buoyancy  flux  B  is  a  key  factor  controlling  the 


Fig.  5.  Separation  of  different  modes  in  the  ( y.p)  plane. 
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modes  of  the  coupled  air-ocean  system.  The  mean 
surface  salinity  flux,  measured  by  y,  makes  the 
upper  ocean  more  buoyant  and  stable.  However, 
the  mean  upward  heat  flux,  measured  by  /a,  makes 
the  upper  ocean  less  stable.  The  larger  the  param¬ 
eter  y  (j a),  the  stronger  the  negative  (positive) 
feedback  mechanism.  Combining  (46)  and  (47) 
leads  to  the  results: 


y  > 


Ip. 


1  -  V/* 


Nonoscillatory  Damping 


2H 


1  -  v/m 
2n>y> 


>  y  >  2/x  =»  Oscillatory  Damping 


2fi 


1  + 


Oscillatory  Growing 


— ^-7^  >  y  >  /a  =>  Nonoscillatory  Growing  (48) 
1  +  yj\i 


Entrainment  regime 

For  the  purpose  of  a  preliminary  sensitivity 
analysis  only  k  and  X  are  allowed  to  vary,  depend¬ 
ing  on  observations.  Other  parameters  are  held 
constant  as  listed  in  Table  1.  By  the  definition  of 
X  (37a)  and  the  estimates  for  3Pr/9n  (9)  and 
Budyko’s  formula  (10),  we  can  estimate  the  value 
of  X  as 

X  -  0.45  X  10-6  s_1 

which  indicates  that  X  has  an  order  of  10-6  s_1. 
Therefore  in  this  paper,  X  varies  from  -10-6  s_l 
to  10-6  s-1.  It  is  reasonable  to  let  k  vary  between 
—  10  to  10.  The  case  of  k  =  1  means  that  salinity  is 
homogeneous  across  the  mixed  layer  base. 

We  compute  all  roots  of  (42)  for  different  val¬ 
ues  of  the  parameters  k  and  X,  and  obtain  three 
roots  at  each  points  of  the  parameter  space  (k,  X). 


TABLE  1 

Standard  values  for  model  parameters. 


p 

1.29kg/m3 

S 

3Sg/kg 

10‘3 

hw 

50m 

hc 

5cm 

a 

0.2  xur 

n 

0.5 

a 

0.8x10* 

*4 

lOm/s 

Cl 

1.0 

Pw 

1035kg/m3 

C2 

0.2 

T 

1  w 

20  °C 

*-h 

1  m/day 

o 


(a)  Re(a,) 


(b)  |Im(ff|)| 


lb) 

Fig.  6.  Distributions  of  eigenvalue  o,  in  the  (k.X)  plane  for 
standard  case  (unit  in  2w/tt):  (a)  Re(o,),  (b)  |lm(a,)|. 


One  root  among  the  three  has  negative  real  parts 
throughout  the  whole  parameter  space,  repre¬ 
senting  the  damping  modes,  in  which  we  are  not 
interested  here.  The  other  two  roots,  a,  and  a2, 
have  positive  real  parts  somewhere  in  the  parame¬ 
ter  space,  representing  the  existence  of  growing 
modes  in  certain  parameter  ranges.  Figure  6  shows 
the  isolines  of  (a)  Re(a,)  and  (b)  |Im(o,)  |  in  the 
(k.X)  plane.  Figure  7  indicates  the  isolines  of  (a) 
Re(a2)  and  (b)  |Im(a2)|  in  the  (k.  X)  plane.  In 
these  figures  the  real  and  imaginary  parts  of  the 
eigenvalues  alt,  a,j,  a2r  an^  o2j  are  scaled  as 
2tt/tt.  Both  the  real  part  (growth  rate)  and  the 
absolute  value  of  the  imaginary  part  (periodicity) 
of  the  roots  a1,  a2  show  saddle-type  distributions. 

Several  interesting  results  from  Figs.  6  and  7 
are  summarized  as  follows: 

(i)  A  necessary  condition  for  the  generation  of 
growing  modes,  which  can  be  seen  from  these 
figures,  is: 

kX  >  0  (49) 
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o 


(a)  Re(a2) 


(b)  1 1  m( a 2)1 


(M 

Fig.  7.  Distributions  of  eigenvalue  o2  in  the  (s.A)  plane  for 
standard  case  (unit  in  2w/tt):  (a)  Refo2),  (b)  |Im(a2)  |. 


This  implies  that  if  the  surface  water  is  warmer 
(cooler)  than  the  deep  water,  i.e.,  k  >  0(  k  <  0),  the 
growing  modes  are  excited.  In  this  case  the  surface 
buoyancy  flux  increases  with  increasing  (decreas¬ 
ing)  cloud  cover. 

(ii)  If  the  damping  modes  are  ignored,  the  two 
eigenvalues  a,  and  a2  exhibit  the  similar  proper¬ 
ties.  The  growing  modes,  a,  and  a2,  are  further 
separated  into  oscillatory  and  nonoscillatory 
modes.  The  condition  for  this  separation  is  given 
approximately  by 

(<  A  Oscillatory 

<cA<  .  . .  ...  (5U) 

(  >  A  Nonoscillatory 

where  A  is  a  positive  number  of  about  2  X  10“ 6 
s  '.  Combining  (49)  with  (50)  and  using  (37a),  the 
necessary  condition  for  oscillatorily  growing  modes 
is 


k  3 
^  >  hc  dn 


PSA  - 


«(^h-  Rs) 


PvT 


pw 


>0  (51) 


and  the  condition  for  the  nonoscillatorily  growing 
modes  is 


k  3 
h,.  dn 


pSsP,  - 


a(Rh~Rs) 

Pw^pw 


>  A 


rw.  n 


(52) 


The  properties  of  thermodynamic  instability  of  the 
coupled  system  depend  largely  on  the  relationship 
between  precipitation  Pr  and  the  cloud  cover  and 
on  the  relationship  between  net  radiation  at  the 
ocean  surface  and  the  cloud  cover. 

(iii)  The  product  of  the  two  parameters  k  and  X 
represents  the  relative  strength  of  positive  to  nega¬ 
tive  feedback.  When  k\  is  larger  than  the  critical 
value  A,  the  positive  feedback  greatly  exceeds  the 
negative  feedback.  The  coupling  system  becomes 
nonoscillatory  and  growing.  When  k\  is  positive 
and  smaller  than  the  critical  value  A,  the  positive 
feedback  slightly  exceeds  the  negative  feedback. 
The  coupled  system  is  oscillatory  and  growing. 
When  kX  is  negative,  the  negative  feedback  ex¬ 
ceeds  the  positive  feedback.  The  coupled  system 
becomes  damped.  Comparing  Fig.  7  with  Fig.  6, 
we  find  that  for  the  oscillatory  growing  modes  the 
two  roots  have  the  same  growth  rate  and 
frequency.  However,  for  the  nonoscillatorily  growr- 
ing  modes  the  growth  rate  relating  to  the  eigen¬ 
value  a2  is  much  larger  than  that  relating  to  a,. 

(iv)  For  the  oscillatory  growing  modes,  the 
growth  rate  or  has  the  order  of  0.5  X  2tt/tt  and 
|  a,  |  has  the  order  of  2-n/rr.  The  corresponding 
period  for  the  oscillation  is 


(53) 


which  means  that  the  period  of  the  unstable  oscil¬ 
lation  has  the  same  order  of  magnitude  as  the  SST 
evoluation  time  scale,  rr,  which  is  20-30  days  for 
the  unstable  atmosphere  and  1-3  years  for  the 
stable  atmosphere. 


Conclusion 


The  feedback  between  the  cloud  and  OPBL  in 
the  coupled  MABL  and  stable  OPBL  system  is 
investigated  by  a  simple  one-dimensional  coupled 
model.  The  time  scales  largely  depend  on  the 
stability  of  the  MABL.  For  the  stable  atmosphere, 
the  two  time  scales  are  quite  long:  rn  T~  100-300 
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days  and  tt  -  1-3  yr.  For  the  unstable  atmo¬ 
sphere,  however,  the  two  time  scales  are  much 
shorter:  t„  t  -  3-6  days  and  rT  ~  20-30  days.  In 
the  western  Pacific  warm  pool  regions,  the  MABL 
is  usually  unstable.  Therefore,  this  theory  may 
provide  some  explanation  of  the  two  time  scales 
(3-6  and  20-30  days)  of  intense  convection  in  the 
western  Pacific. 

In  the  coupled  system,  the  fresh  water  influx  at 
the  ocean  surface  due  to  the  excess  precipitation 
over  evaporation  is  a  damping  factor  (negative 
feedback).  However,  the  surface  cooling  is  a 
forcing  factor  (positive  feedback).  The  relative 
strength  of  these  two  surface  fluxes  determine  the 
mode  type:  decaying  or  growing,  oscillatory  or 
nonoscillatory. 
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